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ABSTRACT 

The sputtering process has been investigated by simulating 
the sputtering of single-crystal copper with 1-7 keV argon. A 
digital computer was used to Bed the crystal, bombard it, and 
move crystal atoms. Four mechanisms Were observed which 
cause surface atoms to sputter. An atom is sputtered when 
(1) it is squeezed out of the surface, (2) it is scooped out 
when another atom strikes its inner hemisphere, (3) it is ejected 
when an atom passes behind it, and (4) it is knocked out by a 
second layer atom which is moving outward. Nearly all sputtered 
atoms were surface atoms. Second and third layer atoms were 
sputtered only for fen energies greater than 5 keV. They were 
sputtered by mechanisms similar to the surface atom mechanisms. 
TSewccwchains" Were OusecrVved=re De Cinected Initem elmer ciel, 
and momentum focuwing was observed to cause sputtering only 
when it occurred in close packed, surface rows. Outward 
directed chains were not observed. Sputtering deposit patterns, 
Sputtering ratios, me sputtered atom energy distributions were 


obtained for (100), (110), and (111) surfaces. All data 
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PN OD Ue Lom 
Material is lost from a substance undergoing ionic bombard- 


Howe. lmh is oemLbemig. Many autm@orc wave eried to define 


Sputtering using terms such as "(atoms)...ejected or knocked 
out"', "emission", "ionic erosion", "disintegration", and “break- 
down." A reader who blithely accepts one of these definitions 


in the course of his reading may conclude: ‘This definition 
describes the mechanism of sputtering. ‘This is an illogical, if 
not erroneous, conclusion; the definition of sputtering must be 
derived from the mechanism, not the mechanism from the 
definition. ‘The sputtering process has been investigated 
utilizing one of two Al mechanisms, ejection or emission. 
Eyjjection implies an meee or direct reaction release of an 
atom from a enyetals emission, a delayed or indirect reaction 
process which results in the release of an atom. 

inwhoe3, Ringdon and Langmuir (1) bombarded thoriated 
tungsten with various ions in a glow discharge tube. This was 
a special case of sputtering since the thin surface film of 
thorium on a tungsten substrate was sputtered rather than the 
tungsten itself. The results of this experiment, that the 
Sputtering ratio, atoms removed per incident ion, increased 
With increased ion mass and increased ion energy, qualitatively 
Suggested an ejection mechanism. A few years later, Von 
Hippel and Eleeneemiac’<? proposed a theory which described 


sputtering aS an evaporation of surface atoms, an emission 


ns 


mechanism. HBarlier, Von Hippel(3) had found by spectroscopic 
means that at least some sputtered atoms were in an excited 
state. The sputtering theory showed that atoms in the region 
of impact could rapidly acquire thermal energy if the kinetic 
energy of incident ions was converted to thermal energy at the 
target surface. If an atom acquired a sufficient area on 
thermal energy, it would then evaporate from the surface, some 
atoms evaporating while in an excited state. This theory was 
improved by Townes(4) in 1944, who calculated a sputtering 
flux and the number of atoms evaporated per incident ion. 
Keywell! 5) in 1954, used neutron diffusion theory to approxi- 
mate atomic Interactions within the crystal, a new approach to 
the theoretical investigation of sputtering. Direct application 
of statistical methods to sputtering was made by Harrison(©) 
wno envisioned the interaction of two distribution functions 
(the crystal lattice and the fon beam). These models, based on 
statistical methods, implicitly accept ejection type mechanisms. 
Onég.of the ~~ important contributions to the study of 
sputtering was made by Wehner!7) im 1953. An the tips 
sentence of his paper, Wehner stated, "The most widely 
accepted sputtering theory Its the evaporation theory. eis but 
the deposit patterns of single crystal sputtering which he 
obtained showed pronounced, high density areas or "spots", 
There was now strong evidence for a momentum transfer 


process, and further development of the evaporation theory 


Ly, 


ceased. Shortly after Wehner's findings were reported, 
Henschke( 8) proposed a theory of sputtering based solely on 
classical collision — treating normal and oblique incidence 
sputtering Separately. Oblique incidence sputtering could be 
explained by two body collision processes. ‘The incident ion 
penetrated the surface layer with little or no interaction, was 
relected outward by second layer atoms, and ejected a sui 
AgFOVMeO yee SCY 1 Lele eelte x its inside hemisphere. His concept of 
normal incidence sputtering required many-body collisions in 
which the ion was Sere reflected outward to sputter 
surface atoms. This theory was plausible for oblique incidence 
sputtering, however, the case for normal incidence LS = uae 
see that the ion be reflected inside the crystal. (For an 
ion more massive than the target this requirement cannot be 
met.) 

Silsbee (9) noted that in any discussion of momentum trans- 
fer effects in sputtering, the geometry of the crystal eho 
be considered; successive collisions in a crystalline structure 
might be influenced by the structure itself. His calculations 
showed: in a two-body collision, an atom has a departure angle 
6; with respect to an axis specified by a close packed direction. 
In subsequent collisions, the departure angle oa will decrease 
if the energy is low enough, and the collision sequence occurs 


along a close packed row of atoms. ‘This concept of "momentum 


focusing", "Silsbee chains", or "focusons! was considered by 


iS 


many investigators, excepting Wehner (10) , to be a Fully satis—- 
factory explanation of the spots in deposit patterns. But the 
deposit pattern is only one of many observable features of 
single crystal sputtering. 

The sputtering ratio is another important characteristic 
of the sputtering process. If a theory or model is to explain 
sputtering, both patterns and sputtering ratios must be 
explained qualitatively and quantitatively. Almén and Bruce 


(AB) ) 


measured sputtering ratios of a variety of metals 
using N, Ne, Ar, Kr, and Xe ions over a 5-65 keV range of 
bombardment energies. (Although polycrystalline specimens 
were sputtered, the qualitative results are the same as those 
obtained for single crystal sputtering.) They noted that the 
sputtering ratio was an increasing function of- bombardment 
energy provided the ion mass was greater than the target mass. 
Sputtering ratios showing a maximum, or of a slowly varying 
nature, were observed in all cases in which the ion was lighter. 
Evidence such as this coincides with any gross conception of an 
ejection mechanism. 

Single crystal Bonner was sputtered by argon at inter- 
mediate energies by Magnuson and Carlson (mc) (12) and South- 


ern, Willis, and Robinson (swr)'13) : 


MC measured sputtering 
Paciosm using W-10 keV argon) for the (111), (100) and»(110) 
surfaces and found that sputtering ratios decreased in this 


order of surfaces. SWR sputtered single crystal copper with 


16 


1-5 keV argon and measured sputtering ratios, but of greater 
significance was the Siena of their deposit patterns wnich 
clearly showed the presence of only certain spots: 

(lit) surface 3-(110), 1-(111) but not seen due to 

beam aperture. 

(110) sw face ee), eer Con 

(100) surface es( 110) ,w be( 100) 

The work of AB, MC, and SWR indicated that the sputter- 
ing ratio is at least a function of bombardment energy (or 
momentum) and of mass ratios. However, consideration of 
spot patterns indicates that the sputtering process itself is 
also highly dependent on the gross and/or surface geometry OT 
the crystal and on ion penetration depth. 

Ion penetration of salide has been investigated experiment- 
ally and by Ben pieermeetniation. Piercy, McCarsco, EPrewi, 


(14) 


and Davies investigating channeling of various heavy ions 
in monocrystalline aluminum, found that the ion penetration 
distance increased in the order (111), (100), and (110) beam 
orientations. ‘This is the same order as the sputtering ratio 
decrease found by MC. A number of studies of ion penetra- 
tion have also been made using computer simulations (Oen, 


(16) 


Holmes, and Robinson‘ 15) , Robinson and Oen Harrison, 
Leeds, and Gay(17)), 
The most apparent correlation between the results of ion 


penetration and sputtering studies is the concept of trans- 


pareney, derined by Elut,” Rieieanel Kistemaker (18), Bott 


157 


experimental and simulation investigations of ion penetration 
indicate that channels exist in certain crystallographic directions 
of a crystal. Results of normal] incidence sputtering experi- 
ments indicates that sputtering ratios are lower for surfaces 
in which these channels are parallel to the beam. One might 
initially suspect that sputtering ratios will be low if ions and 
target atoms are confined to these channels; momentum is 
directed into the crystal rather than laterally, reversal of 
momentum to produce sputtering not occuring. Accordingly, 
some incident ions should see a transparent surface rather 
than a wall of atoms. 

The obvious complexities of the sputtering process may be 
investigated separately by digital computer simulation. This 
of course involves the selection of a gross mechanism (ejection 
or emission) and a model (hardsphere or otherwise). If 
factors such as the interatomic potential function and its 
parameters are known, then the use of the right model in the 
computer may be expected to simulate the actual sputtering 
process. Explanation of deposit pattern characteristics and 
sputtering ratios will necessarily follow from an exact 


Simulation. 


Re 


2. PURPOSE @aanhir INVESTIGATION. 

The belief that the formation of, spots in sputtering deposit 
patterns is Solely a poteeducnee of momentum focusing was 
widely accepted for a number of years. The importance of 
momentum ee oteine. Moony it SPO f OL vert oriymemisicteacmmenien 
entire sputtering process is now believed to be much less than 
originally thought (18-25) , This investigation was undertaken to 
explore the single crystal sputtering process using a digital 
computer to simulate a copper crystal which is being bombarded 
with a beam of argon ions. Primary effort was directed 
towards: 

a. Finding the mechanisms which result in atoms being 
deposited in the spots or high density regions of the 
deposit pattern. 

b. Determining the correlation between crystal structure 
and sputtering ratios and patterns. 


c. Investigating the surface binding energy and its effect 


on deposit patterns. 


ue 


Pome Ko (> Eo... 

A. A space lattice is established whose sites represent 
the equilibrium positions of copper atoms in a face centered 
cule crystal. Atomic spacme is that determined by x-ray 
erystallographic Situdies (for copper, a = 3.015 A). lig y= 
potential function, with parameters for copper-copper ical 


(26). 


action, is the Born-Mayer type Gibson Number ‘T’wo The 


interaction between argon and copper is similarly described 


jen, but the numerical 


using a Born-Mayer potential, V(r) = 
ValiieStOmEne Deeainionows “Suand “peice trhis function are those 
determined by Harrison, een. and Magnuson (HOM) (2) 
from asiguady Gi secondary electron emission. The range of both 
potential functions is eroded at r,, one half the nearest neigh- 


bor ‘distamce sucheihat the potentials and force EO COeZzere=for 


atoms whose centers are farther apart than 2 Use of this 





eroded form of the potential functions permits establishment 
of the crystal in a stable, static state without use of an 
attnraietiine potentials: | This equilibrium state is destroyed only 
when an atom moves from its lattice site. 

The lattice as a whole is restricted by the current computer 
program to perfect form. Defects such as interstitial atoms 
cannot be used, however, irregularities in the surface layer are 
introduced by using two types of variation from a perfect 
SumrTace: vacancy and stub conditions. The perfect or regular 


em@erace is identical to a full lattice plane. In the vacancy 


aw 


configuration the atom which would normally be hit first by 
the — iS ememedetrom the perfect sunrrace, The stub 
surface has an atom added on the perfect surface in a stable 
position and adjacent to the target atom. Additional atoms 
can also be added on the surface or removed from the perfect 
BTS to provide a variety of random irregularities. 

B. Each ion of the beam is approximated by a single, 
neutral argon atom whose velocity vector intersects the plane 
of surface atoms at an impact point. (The term ion is used 
Hho one this thesis for the incident particle to avoid con- 
fusion with target atoms.) 

Successive runs are made using different impact 
points to simulate use of the entire beam area. ‘The use of 
this simple model for the ion beam is subject to four 
conditions: 

1. The results obtained using one random set of impact 
points are essentially the same as those obtained 
using any other random set of impact points. 

a. ‘Vhe majority of beam ions are neutralized before 
reaching some arbitrary boundary which defines the 
surface. 

3. Frior to neutralization, the path of a beam ion as 
not appreciably affected by the surface potential. 

4. A crystal region is in an equilibrium state each time 


a neutralized ion impacts in that region. 


yea 


mach of these conditions has been satisfied either 
in the model itself or by calculation. A full discussion is 


contained in Appendix A. 


rs 


h. DYNAMICS. 

The force functions for copper-copper and argon-copper 
interactions are derived from the respective potential functions. 
Atom position and velocity are then determined using Newton's 
Second Law of Motion. Normally, one would determine these 
values by numerical integration methods over small time inter- 
vals, however, this process uses too much computer time. 
Special methods approximating integration have been developed 
which are used to minimize computer run time yet maintain a 
good approximation to the integration. 

A. Gay(<8) developed an Motaeie method of solution which 
is similar to that used by Gibson, Goland, Milgram, and 
Vineyara!£6). Both Gay's and Gibson's methods use the 
technique of replacing differential quantities with finite 
differences to approximate the equation of motion, F/m = dv/dt 
m/ T! by E/m = Av/eat. Using this equation and the relation- 
ships Ax = vAt and v = (v + vo)/es equations. (3-1) and (3-2) 


are obtained. 


Av=zve-Vv, = | E (x)/m| Noe Ps=1) 

Ee ee = {{F (x,)/m| At/2 + v4 ie (3-2) 
The known values E(x); se x and \t in these equations 
are usually associated fe time t ae PeeaniGe vi, Cle 


unknown values, are to be evaluated using these equations at 
time t = ee (Sieg ae 
Gibson's technique associates the value of XS with a time 


a 


a3 


t =t , but vis associated with a time t = t - At/é. 
ox ° O OV O° 


Cmprent values of .x and .v are alternately computed@a= intervals 
A paising. equations, (3—=™") and (3=2). 


v(t. + At/2) = Mies (to -~At/2) + (F(x, )/m] at 


O° 
x ( tale Ata x(t.) +} [F9)/ ai] at/a + v(t. roveilas 
The advantage of this technique over that of evaluating both 
variables at time t ™ t, + At is that a continuous smoothing 
of the, values of x and occurs with each successive computa-— 
tion. 
Gay recognized the need for Some sort of smoothing 
procedure but considered — an averaging process for equation 
parameters rather than computed results would result in a 


{ 
better approximation. Wis method replaces the evaluated 


force function F(x,) in these two equations with an arbitrary 





force function f(x) which is linear in some interval to be 





determined. Consider the following two step cycle: 


Move an atom from x, to xX, using computed values of 


O° 
F(x,) and vg to solve equations (3-1) and (3-2). Compute 
F(x, ) and average this value with E(x). Now, move the 


atom from x4 to xo using the averaged Ries Cer F(X, x, ) 
and vy, to solve the equations using the same constant value 
of At. Let fnow be defined in the interval (3c x5), such 
iebietosce (eee. ) F( Xi x.) Equations (3-1) and (3-2) may 
now be written: 

vit + At) » v(c) + [F(x5)/m| At 


‘2 Hh. 


mM Coreet) = x( Cc) fel y/aacss 4 wel here 
ie=valics” or yancdex arc determinable in any interval in 
wiieh f approximates = to seme arbitrary degmee, therefore, 
in every complete cycle or timestep a new interval and new 
PUne tiene 1. ware cet tocar 

A determination of which of these two methods is better 
must be based on some standard or condition. In many cases 
Elie aecuracy Of mumerical solutions such as these May be 
compared to the exact value one obtains by integration. ey 
there are no truly exact results in a Simulation since a physical 
model is in the computer rather than the computer being used 
as a means of solving equations. However, if the model is 
physical, physical laws must be satisfied, and the energy gain 
or loss due to the mathematical approximations must be small. 
This condition of energy balance has been used by aeons 
who found that, in simulations involving atomic interactions such 
as channeling or sputtering, use of Gay's method resulted in a 
W-73 energy baiance. 

B. The value of these approximation techniques may be 
substantially reduced by injudicious Selection of a numerical 
value for At. Too large a value invalidates the approximation, 
but too small a value increases computer running time. The 
program used in this simulation of sputtering incorporates not 
only Gay's approximation method but also a procedure which he 
developed for automatic adjustment of At. Its value is 
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calculated using the relationship St = dr/v Where Qt is replaced 
here by bt to indicate that its value varies with each timestep. 
aes velocity,»~y is the magnitude of the velocity of the wmost 
energetic atom in that timestep, and ér is used here as an 
input parameter, the timestep smuleiper. ~This panameee asic 
assigned the computer variable name DTI and is defined: the 
maximum displacement of any atom in any timestep shall be the 
numerical value of DTI in appropriate units of length. Gay 
considered the value of At, consequently that of 1 ee Oud C WE 
function of ion energy, impact point, sand ‘start polne aout 
he was unable to find the relationship. Trial and error methods 
were used to find optimum values, and reasonably good results 
were obtained using these values. Johnson! 9) reported that 
the energy balance was maintained to within 3% for all his runs. 
The difficulties in choosing a value for DTI have been 
largely eliminated in this study. Some problem =r such as 
start point and impact point dependence were avoided by always 
positioning the ion tangent to the first atom !{t would hit. 
This, at worst, will result in a constant error for which 
adjustments may be made (this has not been necessary). The 
positioning procedure is described in Appendix B. The depend- 
ence on ion energy is inherent to the DTI computation process 
and is easily explained once the clutter of impact point: and 
start point dependence has been cleared away. Use of the 
DTI process excludes all variable dependence except: the Inverse 
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proportionality to velocity, but the interval in which the force 
function is considered linear is exactly the same width as the 
numerical value of DTI. Accordingly, a straightforward 
analysis is made to show the interdependence. 

The fom@e fumetion, Er), is expanded about an arbitrary 


point r* using a Taylor's series for a function of one ee 
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Since F(r) is derived from the potential function V(r) = ,AtBr, 
om 1 r 2 “ ° > k 
ss . E(r) is necessarily of exponential form and Fe BE 


1 
Ream and the series is rewritten in the following form: 


F(r)= [even ja + B(e-ee)] Peer (Pr) 2. Se (rone)"| (3-3) 
gas a \ 


The quantity (r-r*) is now defined as a. variable displacement 
Or, and equation (3-3) may be written: F(r) = £(3 r) + g( or) 
where £(0 r) and g( a) are the terms in braces in the force 
function expansion. ‘The function eel ae) is assumed to be the 
linear function ee in Gay's approximation method; the function 
fal é r) is to approximate F(r). Thus, it is required that the 
ratio:\g( § nr) /F( 6 r) be some fraction less than one. 

Let an assumption be made that the ratio is much less 
than one, but that the actual value is the fractional Re cal 
from linearity, a variable whose value is to be specified. ‘Then, 
since B is a negative constant (from the potential function), 


the truncation error for the Series is easily found, and the 
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quantity RO r is known for every specified value of fractional 


deviation. But B is constant, and once the fractional devia- 
tion is specified, dr is fixed. Consequently, specified to the 
fractional deviation from linearity, there is a function flare.) 
which is a valid approximation to the function F(r) in every 
interval (r, r+ dr) where dr is fixed. ‘The quantity dr is 
thus the eine to be used for DTI, and st BC STL LZ ae 
It is now assurred that At will always be the maximum 
possible value which minimizes computer running time without 
introducing errors due to the non-linearity of the true force 
function. One would expect therefore, that the energy 
balance would be maintained to the same order as the linearity 
or@the feree funetion. efurther consideration: will_show that 
the energy balance is maintained to a considerably better 
degree than that anticipated. First, since DTI is the maximum 
displacenfenc™or the most snergetbiesatonm sssentiallywall”otier 
atoms will be displaced a distance less than DTI. ‘Thus, the 
interval in which the force function f(r) is used is smaller, 
and I is aiweven We ttereapproxiniaitienmro” Ber eccond et o:7 & 
given fractional deviation from linearity, the value of §r is 
total displacement, not x, y, or zg component displacement; 
and the velocity used in the determination At = DTI/v for 
elec tiem op isthe tema) velocityey, not ‘the component 
velocity vj. ‘The result of these considerations may be 


summarized by a calculation for the maximum displacement ot 


ae 


am atom along the 1 component of its displacement vector: 


Ax, =v, Qt, butv, Eve and = Say ee Leer ollows 


Vmax 
chai = aE) < DDL. 

These analyses appear to provide a tidy solution to the 
time optimization and energy balance problem. ‘This would be 
true except for the complications introduced by the eroded 
form of the potential function. Consider a situation in which 
one atom is Separated from another by a distance (2r , + qd) 
Wmege ad Is some diatance léSaeithan the total distan@e» D' gehat 
the atom will move in the next timestep. When the atom is 
meved sche tdistamce D,as shown in figure 1, the foreegen tule 
atom at its new position is the same force it would have if 
it had moved only (D-d) units. Since the force is zero for 
py me che medeleassumes chat: Eneratem hadwar Velecin as me 
(D-d)/Ot (which is smaller than its actual velocity v = D/At r 
and a smaller kinetic energy is computed using this smaller 
velocity. ‘This difficulty was originally foreseen by Gay who 
included corrections to the force calculations for such situa- 
wrons. Both Levy! 30) and Johnson! 9) improved the methods 
of corrections, and recent improvements in the force calcula- 
tions have been made by Harmisonee) ° 

The continuous improvement os the dynamics section of 
the program has resulted in a model in which not only has the 
computer running time been minimized, but of greatest Signifi- 


cance, for nearly all of the 700-800 combinations of surfaces, 


EI 


ion energies, and impact points used in this study, the energy 


balance was maintained to within 1%. 
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5. Re aren 

Reports of sputtering studies often include ad hoc formula- 
tions of mechanisms which are used to explain deposit patterns, 
sputtering ratios, and sputtered atom ene distributions, 

The sputtering process is discussed in this thesis in terms of 
mechanisms which have been observed in the simulation to sputter 
atoms. The observance of these mechanisms is an advantage 
peculiar to a simulation. Each crystal atom (and the ion) must 
be identified by number, at least for purposes of computer 
calculations. Since every atom is identified, its complete track 
can be plotted and labeled. The tracks of selected atoms can 
then be superimposed to show a complete set of interactions. 
The mechanisms which were obServed to cause an atom to 
Spueser are conslacred Che prime observable quantity of the 
Simulation. 

The results of the simulation are discussed for each crystal 
Swrace, are Cecer inaction that ail-atoll iaville Citar smers 
sputtered, was made using a probability-of-sputtering PuUNreeron , 
P(E) = 1 - e-E/Ep where E,is the assumed binding energy. 
Wie Vaile OMe ws: compared Wirth a2 random mumber Samad aca 
value between Oandi1. If P)>R, the atom is considered 
sputtered. (This method of selection has been compared with 
one using avstep function in which all atoms with E> eae 
assumed sputtered. A selection method using another probabil- 
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ity-of-sputtering function, P(E) =1-e b? has also been 


tested. No differences Were observed in the sputtering 
patterns based on the atom seleetion method. ) 

Mechanisms observed te eause sputtering are discussed 
first for each surface. Deposit patterns, sputtering ratios, 
and sputtered atom energy distributions are then discussed, 
The deposit patterne shown aré superpositions of these 
obtained by sputtering the erystal separately with regular, 
vacancy, and stub surface conditions. This praetice follows 
that used by Jonneon! #9) whe feund that pattern features 
were difficult to recognize when viewing separate patterns 
from each surface condition, The energy distributions are 
also superpositions of the three surface eonditions. The use 
of superposition is considered to most realistically simulate 
the condition of the erystal surface at various times during 
aputtering. There are strong arguments however, fer using 
only the regular surface eondition, Sputtering ratios are 
not averaged (which would correspond to the superposition 
of other data) sinee a probability faeter for the existence 
of each surface condition would need to be determined, [n= 
stead, the simulation sputtering ratios for vaeaney and stub 
surface conditions are ineluded and diseussed briefly anly for 
completeness. (The (110) surface was sputtered enly at 1 
and 3 keV using aj! three surface conditions. ) 

Crystallographic nomenclature ig used generically in the 
discussion of spots, Referenee te an (hkl) spot dees net 
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imply {hkl) ejection. The¢hkl> directions are with respect to 
a right handed coordinate system in which the crystal is 
described. The (hkl) surface is contained in the x-z plane; 
the y direction is into the drystal. The three crystals used 
in the simulation are shown in figures 2-4; all contain 150 
atoms. ‘This crystal size has been found sufficient to contain 
nearly all energetic collisions for ion Senge ee up to 7 keV. 

7. Whew 11) surface. 

The (111) surface was sputtered normally with argon at 
1-5, 7, 10, 20, and 40 keV. (Although the potential function 
is considered valid only for ion energies less than about 7 keV, 
the higher energy runs were used to search for additional 
mechanisms.) F'our mechanisms were found. These were 
especially evident at the lower (1-5 keV) ion energies. Three 
of these are classified surface mechanisms Since only surface 
atoms were found to participate in the sputtering event. ‘The 
fourth mechanisms is a deep mechanism in which an inward 
moving atom is reflected and then initiates a sputtering event. 
In all cases, the dominant mechanisms were found to be 
surface mechanisms. 

Except for head-on collisions, the ion is scattered by the 
target (atom 2) with a component of momentum parallel to 
the surface as well as a much larger perpendicular component 
into the crystal. ‘The effect of this combination of momenta 


is to drive a nearest neighbor (n. n.), for example atom 6 in 


oS 


figure 5, into an apparent (12> ehannel formed bymaitoms 7, 
8, and 9. ‘The channel is apparent since it terminates 
abruptly thus causing atom 6 to drive atoms 7 and 8 into the 
crystal and atom 9 outward. Since atom 6 passes nearly 
directly behind atom 9, the impulse is more normal to the 
Soe ta cew than pas ela Doom toomedaa moles mechanic 
Since one atom burrows between two layers parallel to the 
surface to sputter an atom in the outermost layer. Atoms 
sputtered by this mechanism are almost always a next nearest 
neigh bomen ant) stousthe target and located in or near the 
sextant defined by the impact area. They usually sputter 
with greater than 10 ev. The target always receives the 
majority of energy transferred by the ion, but it is normally 
driven into the (111) trigonal array of atoms directly behind 
it. These three atoms act as buffers and dissipate the 
Eanoet Sseniewey sd. momentum intomiies crystal. Uhewtanset 
is reflected, but it does not retain sufficient energy to 


sputter a surface atom nor sputter itself. 


The second and third surface mechanisms occur when the 
impact parameter is about one third of an atom radius or 
less. In both mechanisms a n.n., atom 6, is struck by the 
ion such that the parallel and perpendicular momentum compon- 
ent magnitudes which it acquires are of the same order. If 
the perpendicular component is greater, atom 8 scoops a n.n.n. 


atom 9, relatively high up in the surface plane of atoms and 
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into their outer hemispheres. These atoms are densely packed 
in Ghe surface and acteas a rigid reflector for atom 9. It is 
DOMiGmetaieacoms Sputtered by wthe Sscocpamechaniaom lave 
ejection angles nearer the normal than would be found in the 
absence of reflection in the surface. If the perpendicular 
momentum component of atom 6 is nearly the same or less 
Eiamemehe parallel component, it squeezes atom 9vagainee Ics 
surface neighbors. ‘The squeezing causes the surface plane to 
warp and atom 9, 10, or 11 is sputtered. Atoms sputtered 
byeeorenicr me Sseccop- or Sduceze mechanism usually "spuecer wien 
less than 10 ev. Atoms sputtered by combinations of the 
three surface mechanisms are found to sputter with higher 
energies than when sputtered by a single mechaateem 

ie Tomren mochanism™ requires that anwavom™ be remocemad. 
The atom which is most often apt to be reflected is the 
target. If the impact parameter is greater than one half 
the atom radius, the target is driven towards an edge of the 
trigonal array behind it rather than into it. It will penetrate 
this edge if it has sufficient energy, but it will lose most of the 
energy during the penetration. ‘The reflection occurs from 
third layer atoms and the target atom's reversed momentum is 
transferred to a second layer atom. ‘The second layer atom 
will then sputter one or two surface atoms at energies upwards 
of 5 ev. At lower ion energies the target atom cannot 


penetrate the array edge and is reflected from the second 
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Iegcore wit omay or may not Sputter a°snretecera tom, see Ic 
rarely sputters itself. 

When the crystal was sputtered at 7 keV, the ion was 
found to penetrate well into the second layer. ‘This is a 
significant increase in ion penetration depth since 5 keV ions 
are found to penetrate only to the order of half an atom radius. 
The effect of deeper ion penetration is to cause a lateral 
compression of the second layer with accompanying warping. 
Second layer atoms which are squeezed outward will sputter a 
surface atom but are rarely found to sputter themselves. At 
the higher bombardment energies tested (10, 20, 40 keV), the 
number of Second layer atoms which sputtered increased nearly 
proportionally to the increase in ion energy. The sputtering 
process at theSe higher energies appears to be predominantly 
by the three surface mechanisms described for lower energy 
Sputtering. Mechanisms peculiar to high energy sputtering have 
not been observed although this may be due to the limited 
energy range for model validity. 

All (111) simulation deposit patters showed the character= 
istic features of three (110) spots and a central (111) spot. 
The 2 keV pattern definitely showed the presence of three 
additional spots which were found in the regions in which 
streaking occurred at other energies. These patterns are 
shown in figures 6-21. The appearance of a hexagonal pattern 


at « keV was a surprising result, but it is not a unique 
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occurrence in sputtering studies. Anderson and enon to) 


found a hexagonal pattern for (111) copper sputtered by 
mercury. The expected trigonal pattern was found for 
energies up to 400 ev, a hexagonal pattern appeared at 400 
ev, and it disappeared as the ion energy increased. Very 
recently, Robinson and Southern! 32) have found additional 
spots near < 114 positions for (111) gold sputtered with 

4h keV argon. 

The pattern features vary slowly with binding energy 
provided its value is restricted to the range 2.50-3.50 ev. 
The upper limit is the sublimation energy for aegis (33); 
ayove Tc tne pattergs appedr to deteriorate. Bilis wacemicse 
evident at a 3 keV bombardment energy. A set of patterns 
at this bombardment energy is shown in figures eee for a 
binding energy range of 1.50-5.50 ev. The pattern deteriora-— 
tion is most evident in the region of the (111) central spot. 
Figures 11 and 12 show the transition from a well defined 
central spot at 3.00 ev to complete deterioration of the 
central spot at 4.00 ev. 

The half-intensity width of the (110) spots is estimated 
from pattern data and from numerical data to be about 11 
degrees. The width of the central spot is an unreliable 
datum and was not determined. 

Sputtered atoms were always found to be surface layer 
atoms for sputtering at 5 keV or less. The atoms which 
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were sputtered most frequently appeared to naturally group by 
energies into two categories; atoms with perpendicular energy 
greater than 3 but less than 10 ev and atoms with perpendicular 
energy greater than 10 ev. The higher energy atoms were 
found to be equally distributed between the three (110) spots 
and the central spot. ‘They were also found in the narrow 
sectors defining the streaks and at distances from pattern 
center corresponding to the distances of the (110) spots. 

This distance was found to be 0.7 units, the distance measured 
for the spots in the patterns of SWR. The lower energy 
atoms were usually found in the central spot which explained 

its Sensitivity to binding energy. 

The correlation of sputtering frequency and crystal 
Ipeaition ofmsputtered atomauni2s giopnd 1Slncet Tr eaiienemalocairion 
diagrams. These are shown in figures 22-29 Poi regular 
surface sputtering. Similar diagrams have — made for 
sputtering the vacancy and stub surfaces. It was found that 
a 1 keV ion does not cause frequent sputtering of n.n.ni's or 
nen.'s in the sextant defined by the impact area. However, 
aiuwe and 3 keV, nen.n.'s are Pound to be sputtered most 
frequently, and at higher energies, both — amidase Ss 
are frequently sputtered. The three-fold relationship between 
the frequency of sputtering, location of the atom with respect 
to the target, and the ion energy may be clarified by consider- 


ing the sputtering mechanisms. 
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At 1 keV and an impact parameter of about three fourths 
an atom radius, the ion is found to penetrate less than half 
an atom radius. At energies up to 5 keV the penetration is 
only slightly greater. It was found that, although the target 
is always driven into the crystal, the n.n. which initiates the 
majority of sputtering may not be driven far enough into the 
Gmetam ve provide the scoop mechanism. If Weiser driven 
ide ametiticicnt distance to scoop the n.n.n«w«, ibsmenersy 
will be propagated in the surface along close packed rows 
which orfginate at the n.n.n. ‘This is not a "focuson."' At an 
early point in the propagation the surface will have warped 
sufficiently to cause an atom to sputter. ‘The sputtering is 
a result of a squeeze mechanism but not one directly involving 
the n.n. As the ion energy is increased the scoop and “mele 
mechanisms are more apt to occur although the squeeze 
mechanism is found to occur at all ion energies. ‘The transi- 
tions from dominance of one mechanism to another with 
increasing ion energies may be inferred from the profiles shown 
in figure 31 for atoms which sputter most CRAG ol be 

A relatively constant sputtering frequency is observed for 
atoms 102 and 132, but the percentage of high and low energy 
Sputs varies for atom 102 °whiech is a°ntnei {hem pereomeag- 
degesmiOt Vaigy for atom 132 which is neither a n.nvne nord 
nen. Atom 102 sputters most frequently at low ion energies 


by. the squeeze mechanisms, at higher ion energies by the scoop 
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and mole mechanisms. These relationships have been €ontirmed 
by the atom track displays. They also have shown that atom 
132 is sputtered by atom 102 through a Squeeze mechanism 
alone at low ion energies but in conjumetion with jehe mole 
mechanism at higher ion energies. The profile shown for atom 
86 is the result of simple reflection from gece layer atoms. 
Atoms 87 and 101, which are not in the impact Sfeporelann loyohe 
are Spubeevicd sawedient)y eae opiueoRed by ache woolce ac 
mechanism. 

Pi Ces acon SOR Omiles Of satan s iS plRiemeda by emciece 
mechanism. At low ion energies the target reflects at an 
oblique angle from Second layer atoms rather than penetrating 
the layer, and it enters one of the apparent ¢112» channels. 
Atoms 25 and/or 55 are then sputtered directly by the target. 
(Although this is similar to a mole mechanism it is a distinct 
mechanism since the target must be reflected in this case but 
not in the former case.) At higher ion energies the target 
will penetrate the second layer and is reflected from third 
layer atoms. It does not channel, but it causes a second 
layer atom to sputter surface atoms. Accordingly, atom 55 is 
sputtered by sthiegdeae Bnecianienmimonme sol ten thansatom 25 is 
Sputtered by a channeled target. Atom 40 or 70 is 
sputtered by a Squeeze mechanism, usually in conjunction with 
the sputtering of atoms 25 and 55. 


These results for the sputtering of a regular surface 
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have been found to be generally applicable to sputtering vacancy 
and stub surfaces. Equal numbers of atoms from the three 
surfaces are found in the spots and in the Streak regions. 

It is also found that these atoms are usually the same from 
all three surface conditions. They differ only in energy since 
it is the ion penetration distance which indirectly determines 
the magnitude and direction of momentum transferred to 
surface atoms. The vacancy surface Sees an ion which 
penetrates deeper; this is analogous to an ion of greater energy. 
In the stub condition, the ion does not penetrate the surface; 
it transfers energy to the stub atom. ‘The stub then assumes 
the role of the ion impinging on the surface, but the stub 
appears as an ion with lesSs energy. 

Pigure 33 shows the sputtering ratios for each surface 
condition as a function of ion energy; assumed binding energy 
is a parameter. A statistical variation of +204 of the 
sputtering ratio has been assumed. ‘This would normally be an 
unreasonably large deviation for laboratory results, but it is 
considered conservative for the numerical results from a 
simulation in which many parameters are unknown. ‘The actual 
value of sputtering ratios obtained for the regular surface is 
generally low by 2-3 atoms/ion for a binding energy of 3.5 ev. 
Further investigation Showed that, in a few cases, surface 
atoms 3 planes below the impact point would be sputtered with 


about h-5 ev. The average increase in the (111) sputtering 
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ratio over the i-7 keV range is estimated to be no greater 
than 2.0 atoms/ion for a binding energy range of 1.5 to 3.5 
ev. This range is shown by the dotted lines in figure 32a. 
Vacancy and stub sputtering ratios show erratic behavior over 
the energy range studied. 

Energy distributions of sputtered atoms are shown in 
figures 35-39. ‘The preponderance of atoms with energies less 
than 5 ev is caused by the lack of an intrinsic surface binding 
energy. The sputtering selection process considers all atoms 
having EB Pe eeiieG than 1 ev. Thus the P(0.50) energy, using 
an assumed binding energy of 3.5 ev, is @.4e ev. This is 
especially evident for the 1 and 5 keV distributions. If this 
region is disregarded, the peaks which occur at 5-8 ev for 
all ion energies are assumed to be the maximum for each 
distribution. Secondary maxima appear at 16 ev for 1 keV 
Sputtering, 16 ev for 2 keV, and at 1h ev for 3 keV sputter- 
ing. At higher ion energies, a secondary maximum may be 
present at 47 - for 5 keV sputtering. The number of atoms 
sputtered in the simulation is too small to make a definite 
statement regarding maxima. Differences are usually measured 
baone,..at Most, <Ewes atoms: 

B. The (100) Surface. 

The (100) surface was sputtered normally with 1, 3, 5, 
and 7 keV argon. Mechanisms observed to cause sputtering 
were found to be nearly identical to the surface mechanisms 
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discussed for (111) sputtering. A mechanism directly 
cmp tO tie deep Mechalisn Was Mem ouGem=d. IMSstead, 
fieewes would that the Scoop and nole mechameameaeare emlanced 
by the presence of the (110) channels parallel to the surface. 
Surface atoms which are driven into these channels do not 
have to burrow between the first and second layer in order 

to sputter nearby surface atoms. Additionally, it was found 
that these atoms will frequently sputter themselves by 
reflection from second layer atoms. ‘The squeeze mechanism 
was observed to be the most effective sputtering mechanism 
at 3 and 5 keV. At 1 and 7 keV ion energies, variations of the 
scoop and mole mechanisms were dominant. The low-high and 
middle ion energy dependence observed for the dominant 
mechanism is directly related to the shape of the impulse by 
which the tion transfers energy to a n.n. 

At low ion energies, the impulse is sufficiently broad that 
the n.n. is directed into the edge of the Square of atoms 
penind it. The n.n. will enter the (110) channel without 
immediately sputtering a surface atom. When the impulse is 
narrow, the peak force is not necessarily greater. The target 


has received most of the energy given up by the ion, and the 


ion—={n.Tl. impact parameter Is larger.” 7 hewiie ie Wile eire emer 
or Squeeze a surface atom causing it to sputter. At high ion 


energies (7 keV) the ion was found to penetrate at least into 
the second layer where it was deflected towards but not into 
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a (110) channel. This was observed to cause second layer atoms 
to be scooped up and sputter surface atoms. Second layer atoms 
are also squeezed by their neighbors and may sputter through 
the vacancy left by a sputtered surface atom. 

Sputtering at 3 and 5 keV was found to be caused by the 
scattering of surface atoms along the surface as the ion was 
being reflected from second layer atoms. The sq ueeze 
mechanism was observed to sequentially sputter atoms along 
a close packed surface row. The occurrence of a mole 
mechanism was conspicuously rare; even at 5 keV the energy 
propagation was clearly restricted to propagation parallel to 
the surface layer of atoms and to propagation Into the crystal}. 
The only definite momentum reversals observed were for the 
target atom and the ion. 

The characteristic features of (100) sputtering deposit 
patterns are four (110) spots and a central (100) spot. The 
patterns are usually outlined by a hypocyclold-shaped haze} 
the (110) spots form the cusps of the hypocycloid. Simula- 
tion patterns were found to show these characteristic 
features quite well for 3 and 5 keV sputtering. The 3 keV 
pattern (figure 41) is a very good likeness to one at: 2.5 keV 
reported by SWR. The (110) spot distances from pattern 
center were found to be 1.0 TC unlt whieh corresponds to 
the distances in the 2.5 keV experimental pattern. It was 
further determined from numerical data that the (110) spots 
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ss actually elliptical rather than circular. ‘This latter 
characteristic is undoubtedly the result of distortion inherent 
in flat plate collection of sputtered atoms (13) | The 3 keV 
pattern was found to remain essentially unchanged as the 
binding energy was varied from 2.50 to 3.50 ev. The 5 keV 
pattern (figure 42) was almost identical to the one at 3 keV 
only when a binding energy near 2.00 ev was used. As the 
binding energy was increased, the (110) spots became less 
well defined (figure 43) and spots corresponding to 2210) 
became the most prominent feature. The intensification of 
¢210) regions was accompanied by a loss in definition of the 
hypocycloid outline. Similar results were obtained when sputter- 
ing at 7 keV. The hypocycloid outline in the pattern could 
only be observed for a binding energy near 1.50 ev. 

The deposit pattern for 1 keV sputtering (figure 40) 
showed very little similarity to those obtained from higher 
energy sputtering. F'our spots corresponding to ¢211> were 
observed; four (110) spots may also be defined but they are 
extremely diffuse. (The (110) spots are best seen in the 
point plot in figure 40a.) ‘The finding of (211) spots at 
1 keV, and the appearance of (210) spots for 5 and 7 keV 
ion energies suggested that the pattern might be rotated 45 
degrees for low ion energy sputtering. ‘The possibility that a 
Peracron Occurs Nas been Mvestigarea (25m but one has never 


been experimentally observed. 
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The different (100) pattern features which appeared with 
varying ion energies and binding energies were a mamled eanieras t 
to the relative constancy observed for (111) patterns. (The 
use of a binding energy of 3.00 ev resulted in comparable (111) 
patterns for all ion energies.) The lack of constancy for the 
pattern simulations may be explained as a result of the sub- 
traction of a binding energy from the Set eran ianliaes component 
of a sputtered atom'!s energy. (This is done to simulate the 
energy lost in overcoming the surface potential. ) es 2m 
sputtered with small perpendicular energy, less than, 15.eyetor 
example, may suffer an apparently small change in its perpendic-— 
ular velocity, but the direction of its velocity may be consider- 
ably waltered. Thesavmevageupeupendiculenwcuicrcy ofa Gomsuheund 
in spot regions is about 15-20 ev from sputtering the (100) 
facees Lisic abou. 8=——>.c cummed ( is) sputtering fonsanadeh 
variations of up to 1 ev binding energy have had little effect on 
pattern features. ‘This explanation does not, however, explain 
why the pattern at 3 keV is valid for a range of binding 
energies. A consideration of the atoms sputtered may clarify 
but not completely explain this situation. The dissimilarities 
seen in the (100) patterns may be directly inferred from the 
frequency-location diagrams for (100) sputtering (figures 6-49). 

When the energy propagation is restricted to the surface, 
as it is for 1 keV sputtering, the (111) surface can dissipate 


energy through four atoms in the momentum forward semi- 
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circle. ‘The (100) has only three atoms comparably located 


to dissipate the energy. ‘Track displays have shown that for 

1 keV sputtering, atoms 26, 116, and 132 are sputtered by a 
combination of mole and Squeeze mechanisms. The colemts 
predominant in sputtering the atoms nearest the target, the 
Squeeze mechanism becoming predominant for surface atoms 
which are located further from the target. Atom 26 sputters 
less frequently than atom 116 since it is on the see) Side oT (“10 
of the crystal with respect to the impact area. Both of 
these atoms are found in (211) spots, but atom 116 is also 
found in (110) regions of the pattern. Atom 132 is found in 
(211) spots and the (100) spot; it is TKS only by the 
squeeze mechanisms. When the bombardment energy is 3 keV, 
the deposit pattern has its expected characteristics. Correla- 
tion of high atom density regions with frequency-location data 
showed that the (110) spots are the result of sputtering of 
n.nen.'s. Concurrently, the hypocycloid outline was found to 
be formed by: (1) the sputtering of the same atoms which 
formed (211) spots at 1 keV, and (2) the sputtering of the 
n.en.n.'s in the surface. These atoms were not densely 
deposited in the (211) regions but formed diffuse ellipses whose 
semi-minor axes were along quadrant biscenenes The appearance 
of well defined (110) spots at 3 keV is strongly suggested by 
the 3 keV frequency-location diagram, figure 47. Atoms 70 


and 72 were observed to be sputtered with about the same 
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frequency as atoms 55, 56 and 86. Since atems 70 and 72 
are along a (100) axis, one might expect: that they would be 
found in (110) spots, and this has been observed. Atoms 55 
and 56 are in positions relative to the target such that they 
would sputter within adjacent 45 degree sectors in quadrants 
TI and [f!f. ‘These atoms are found In the hypocyeloid outline. 
The absence of the hypocyeleld at 1 keV and lts presence 
at 3, 5 and 7 keV Is readily apparent when one eonsidei's the 
sputtering profiles in figure 50. At i keV ton energy, n.n.'s 
are never sputtered, however, they are sputtered with 
increasing frequency 48 the ion energy !nereasges. Atom 86, 
the n.n. In the quadrant containing the impact areu, exhibits 
this behavior of increasing sputtering fFrequeney., Atom 72 
shows 4 sputtering profile which indleates the (110) spot 
formation at lon energies of 3 keV or greater. The constant 
nature of the profile tor atom 102 at 3 keV and higher ton 
energies is indicative of the hypowyelold outline rather than 
C211) spot formation. Atoms 26, 116 and 132 show 4 generally 
decreasing frequency of sputtering with high ion energy. ‘This 
is consistent with the loss ef (211) spets at higher ton energies. 
The track patterns observed fer 5 and 7 keV (100) 
sputtering showed that atoms ejeeted by 4&4 Squeeze theehaniam 
are often ejected in directions appesite to those ane would 
expect solely on the basis of the loeation of the atam with 
respect to the Impact point. Tiedr eeamr baer Were @atised 


1, & 


by two factors: (1) the atom was squeezed against its 
neighbor and reflected outward and with its parallel component 
of momentum reversed rather than being strictly ejected out— 
ward. (2) Sputtering by near simultaneous squeeze and mole 
mechanisms results in ejection of an atom where the ejection 
direction is dependent only on the impulse delivered by the 
channeling atom in the mole mechanism. 

Sputtering ratios for the (100) surfaceare shoWmin figure 51. 
The correspondence between simulation values and experimental 
values is quite good for the regular surface for a binding 
energy of 3.50 ev. The ratios for vacancy and stub surfaces 
show a closer correlation to the regular surface for (100) 
Sputtering then tor (Ti1)"sputterine. 

The energy distributions of sputtered atoms are shown 
in figures 52-55. No specific maxima other than that near 
5-7 ev appeared to be present. 

Cr net 11) “Surrace. 

This was the last of the three face centered cubic crystal 
surfaces to be sputtered in this simulation. It is purposely 
the last to be discussed. First, experimental deposit patterns 
from (110) sputtering show only a large central oval area for 
bombardment energies greater than a few hundred ev (10). 
Second, a unique feature of this surface is the (110) channel; 
no other (hkl) surface in a face centered cubic crystal shows 


Chkl> channels. The sputtering mechanisms OCCUrring@rere sy 3, 
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5 and 7 keV ion energies were most readily observed iene 
sputtering of this surface, and they are interrelated —_— the 
presence of (110) channels. The mechanisms are, again, 
identical in concept to those previously discussed for (111) 
and (100) sputtering. Before discussing these mechanisms 
separately in terms of individual crystal atoms, the features 
of the frequency-location diagrams (figures 56-59) are 
summarized. 

Sputtering of the (110) surface differs significantly from 
the sputtering of (100) and (111) surfaces: (1) the target 
atom was the atom sputtered most frequently at al) ion 
energies. (This is in agreement with results obtained by 


(30) ) (2) the number of sputtered atoms with high 


Levy 
energies were generally those in the close packed row contain- 
ing the target atom, but with the. target atom located at 

the center of the row rather than at the origin of the row. 
These findings are wholly consistent with expected results 
when one considers the sputtering mechanisms with respect to 
this particular surface. 

The trequent sputtering otf the target is made possible by 
the nature of the squeeze mechanism but the mechanism does 
not itself cause the sputtering. When the target is driven 
into the crystal it squeezes the n.n. (in its row) and ejects 
it. Once the n.n. position is vacant, the target will travel 
through a large, potential-free area betore striking a second 
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layer atom. ‘The target strikes this second layer atom and is 
free to reflect outward without finding a surface atom directly 
in its path. The target then transfers most of its parallel 
momentum to an.n.n. in the row. The parallel impulse is 
then propaéagated down the surface row. Warping of this one 
row causes additional atoms to sputter. ‘This sputtering 
Sequence ts also observed as a result of the lon sStrikimeecie 
n.n. which is in the quadrant defined by the impact area. A 
rather surprising continuation of this sputtering mechanism 
was observed to occur in atom rows located both above and 
below the horizontal row containing the target. Second layer 
atoms which receive energy from the target or ion are driven 
into (111) apparent channels. They pass behind another second 
layer atom causing it to be ejected outward and strike two 
adjacent surface atoms. Sputtering is, again, initiated in a 
surface row. Sputtering profiles for the target (atom 72) 
and atoms h2, 71, 73 and 102 are shown in figure 60. Atoms 
42 and 102 are the n.n.'s to the target; atoms 71 and 73 are 
those which initiate sputtering in the horizontal rows containing 
them. Second layer atoms are sputtered more frequently as 
the ion energy is increased. 

Simulation deposit patterns from (110) sputtering (figures 
61-64) tend to show more of an oval outline than a uniformly 


dense central oval region. The central region can be made 


more dense by including atoms sputtered with less than 1 ev 


oy: 


perpendicular energy in the probability-of-sputtering selection 
process. This did not seem justifiable for (110) pattern 
production since these small energy atoms were excluded Irsrouge 
(100) and (111) patterns. The 1 keV point plot (figure 61a) 
showed indications of (100) spots but these cannot be clearly 
seen in the deposit pattern (figure 61b). Their presence was 
substantiated by numerical data. The average energy deposit 
data showed that the average energy of atoms in 1005 regions 
was nearly double that of atoms found in the central oval. 
This is a constant characteristic of non-central spot regions. 
At higher ion energies there was no direct indication of the 
presence of these spots either by numerical data or energy 
deposit data. A few atoms having high energy were found in 
these regions but the area density was too small to form a 
spot in the pattern. The set of atoms found to form the 
central oval for sputtering at all ion energies was not a well 
defined group such as has been found for (100) and (111) 
Sputtering. The atoms forming the (100) spots trom i keV 
sputtering are those which are either n.n.'s to the target 
or mvun.'setom@thesaitom which initiatzesmeputtering in ceach 
Momizontal rows 

The sputtering ratio for (110) sputtering is shown in 
figure 65. Agreement with experimental data from the 
numerical and curve shape aspect: was poorer for this surtace 
than the other two surtaces. ‘his is apparently caused by 
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(110) sensitivity of the sputtering ratio to the A parameter 
(30) 


in the argon-copper potential function. Levy obtained 
better curve shapes using A = 11.435 whereas A = 12.56 
has been used in this study. This sensitivity has not been 
observed for (100) and (111) sputtering. 

Energy distributions are shown in figures 66-69. The 
notable feature of (110) sputtering is that the energy of 
sputtered atoms is generally higher for (110) sputtering than 
(100)or (111) sputtering. 

D. Results of Ancillary Studies. 

During the course of the simulation, it became apparent 
that certain extensions of this study should be made. It was 
particularly desirable to probe five areas: 

(1) Sputtering mechanisms peculiar to high bombard- 
ment energies (up to 40 keV). 

(2) Sputtering of the (110) surface by high energy 
ions which are directed towards (110) channels. 

(3) Sputtering mechanisms for an ion heavier than 
the target. 

(4) Sputtering of the (0001) basal surface of a 
hexagonal close packed crystal. 

(5) Sputtering the (100) surface at ev ion energies. 

Investigation of the first two of these five areas were 


limited somewhat by the restrictions for the copper=—copper 


and argon-copper potential function. ‘They are usually assumed 
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valid up to about 10 keV, however, if small impact parameters 
are assumed to rarely occur, valid results may be expected. 
Areas (3) and (4) involve the use of unknown parameters. 
Accordingly, the validity of results from investigating these 
two areas cannot be assured. 

The search for mechanisms peculiar to high energy sputtering 
was made using the (111) surface. This surface was considered 
the one most likely to show additional mechanisms since it is 
the most densely packed surface. ‘The failure to find mechanisms 
(part A of this section) peculiar to high energy sputtering 
was not too surprising. 

The sputtering of the (110) surface by ions directed into 
(110) channels was investigated by determining the sputtering 
ratio for ions impacting at the two impact points located 
nearest the channels. The sputtering ratio curve for channel 
shots is shown in figure 70. ‘The discontinuity between the 
low and high ion — curves is assumed to be the result of 
breakdown of the potential functions. The important result 
is that sputtering always occurred for channel shots, even 
when 40 keV ions were used. 

The third area was examined using xenon to sputter the 
(111) surface. Argon-copper potential parameters were used 
Since those tor xenon-copper were not available. Runs were 
made at 3 and 5 keV. At both energies, the ion penetrated 
agelesaSeebemthe fourth layemeof atoms;,no,indicationof momentum 
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reversal of the ion was obServed. ‘The ion was found to be 
channeled, just before penetrating the third layer, into (110) 
Ee cite It initiated cascades which propagated into the 
crystal. 

The fourth area was investigated to see if the effect of 
surface geometry on the sputtering deposit patterns could be 
determined by a comparison of (111) sputtering of a face 
centered cubic crystal with sputtering the (0001) surface Or 
a hexagonal close-packed crystal. Hasiguti, Hanada, and 
Y cio se? have sputtered zinc with 8 keV argon; the 
deposit pattern showed an outlined, equilateral hexagon Ww eine: 
a central haze. An attempt was made to reproduce this 
pattern with the simulation model modified for zine~zinc 
interactions. The potential parameter was adjusted for a 
Born-Mayer type potential V(r) = Aer? Where A = 52 
(ayz,, 0h keV. ‘This relationship was determined by Ander- 
sen and Sigmund! 37). The argon~-zinc potential parameters 
(assuming a Born-Mayer type potential) were approximated by 
using those for argon—copper interactions. Neither the 
sputtering ratio nor deposit pattern matched the reported 
results. This was not unexpected since the zinc atoms in 
the crystal must be represented as ellipsiods (with the minor 
axis in the basal plane) rather than represented as spheres. 

The (100) surface was sputtered at 100 ev to investigate 
the possibility of a 45 degree pattern rotation tor loweaiem 
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energy sputtering. ‘The number of atoms which were sputtered 
was extremely smail. Ali sputtered atoms had energies less 
than 10 ev. It was found that the atoms sputtered most 
frequently from the regular surface were n.n.n.'s (figure 71). 
Bombarding vacancy and stub surface conditions also resulted 

in a predominance of n.n.n.'s being sputtered rather than 
sputtering of the atoms forming the square about the target 
atom. 

E. Results Common to All Surfaces. 

Momentum focusing was observed to oceur onl, integthe 
crystal for the three copper surfaces studied. This was 
particularly evident from atom track displays of (111) surface 
sputtering. If an atom underwent momentum reversal, it was 
always by reflection from atoms located no deeper in the 
crystal than the third layer. Ewen for 7 keV ion energies, 
crystal atoms located deep in the crystal were always driven 
inward. A Surprisingly large number of atoms with large 
energies were found to move between second and third layers 
with their motion nearly parallel to the surtace and with a 
small, inward-directed momentum component. ‘lhe energy 


which appeared to be delivered in impulses to atoms still at 


lattice sites, dissipated through the crystal. ‘These atoms 
rarely acquired more than 20 ev through the energy impulse 
process. ‘l‘heir motion was restricted since they were 


surrounded by other atoms. 


The target atom in the (100) and the (110) surfaces, 
when struck at near zero impact parameter by the ion, 
always transmitted the majority of its energy billiard—ball 
fashion in close packed rows perpendicular to the surface. 
(Atoms to which this energy was transferred escaped through 
the back face of the crystal.) ‘There were a few cases in 
which one of these atoms would initiate another chain. ‘The 
chain was never more than 2-3 atoms in length before the 
energy was either wholly dissipated or the chain stopped by 
divergence of the momentum to form numerous small cascades 
directed into the crystal interior. 

The unique arrangement of atoms in the first layer was 
observed to be a dominant factor in determining ejection 
directions for all surfaces. Atoms which were ejected other 
than nearly normal to the surface were always influenced by 
their neighbors. Atoms which were ejected at angles near 
45 degrees to the normal would reflect from their neighbors, 
often ending up as normally ejected atoms. ‘This effect was 
most pronounced for (111) sputtering since the six atoms 
acted as a lens, but it was also seen for (ire) sputtering 
in which a lens is formed by second layer atoms. The (100) 
surface showed a strong lens effect although one would not 
necessarily assume that it occurs. Sputtered atoms, 
originally adjacent in the (100) surface, were often observed 


to be ejected almost simultaneously. In these situations they 
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were ejected nearly normal to the surface. This may be the 
genesis of dimers recently observed by Woodyara!32) , 


The spot regions of the deposit patterns from (100) and 


(111) sputtering always contained more high energy (> 10 ev) 


than low energy atoms. But the energy distribution was more 
uniform than increasing with the ion energy. It was not 
peaked at any one or group of energies. Atoms with low ion 


energies were predominantly found in the central region of the 


Dawe mii. 
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6. Conclusions. 

The observable quantities of sputtering appear to be 
interdependent on only the sputtering mechanisms. There is 
no evidence that these quantities are interrelated. The 
qualitative and quantitative data from the simulation indicate 
that the deposit pattern, sputtering ratio, and energy 
distribution of sputtered atoms cannot be correlated with 
each other; one cannot predict a sputtering ratio from an 
energy distribution. But, each of these quantities can be 
cross-—correlated between surfaces. 

The main features of deposit patterns appear to be 
determined only by the surface geometry of the crystal. ‘The 
formation of (111) pattern spots is attributed primarily to 
assisted focusing by the hexagonal lens. The predominance of 
trigonal rather than hexagonal symmetry is considered a 
natural result of the brief (112) channeling observed in the 
mole mechanism. The appearance of streaks between the 
Spot pairs and the appearance of a hexagonal spot pattern at 
2 keV suggests that use of a hemispherical collector in the 
simulation will show the three (114) spots which were observed 
by Southern and —epemen coe Fiowever, no distinction can 
be made between the (110) and (114) spots on the basis of 
pattern location. If the mole mechanism is predominamigmae 
certain energies the (110) spots should be more intense than 


the (114) spots. If the squeeze mechanism is predominant, 
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the (110) and (114) spots should have more equal intensities. 

The nearly exact simulation of a (100) pattern at 3 keV 
is considered one of the best arguments that sputtering is 
mostly a surface phenomenon. Neither the potential form 
nor its parameters for copper-copper and argon-copper inter- 
actions are known with certainty. If sputtering was a deep 
phenomenon, one might accept an argument that spots in the 
pattern could be produced without exact knowledge of the 
potential. However, the hypocycloid outline, definitely present: 
in the simulation pattern, is considered the feature of the 
pattern which would be most sensitive to small variations in 
the potential if sputtering involved more than the first few 
layers. ‘This conclusion is consistent with the apparent 
rotation of the 1 keV pattern. Atoms found in Ci regions 
were sputtered when the energy transferred by the ion to the 
target was as small as 40 ev. This compares favorably with 
a calculated transfer value of 75 ev for a 0.5 A Ween crsnene 
parameter (21), A sputtering threshold energy of 50 ev 
has been reported (32), 

The inability to produce consistent (100) patterns at a 


given binding energy for 5 and 7 keV sputtering is undoubtedly 


due to the unknown condition of the surface. It would be 
unreasonable to assume that the binding energy is a decreasing 
function of ion energy only for the (100) surface. ‘I'he 


disruption of the surface by the first group of incident ions 
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is a factor which cannot be ignored. ‘The sputtering ratios 
determined for the (100) surface as well as for the (110) 
and (111) surfaces, were within reasonable limits when the 
assumed binding energy is in the range 1.50 to 3.50 ev. No 
definite value of binding energy can be determined from the 
patterns and sputtering ratios unless a weighting factor is 
used for the surface condition at the time of ion impact. 
The assignment of weighting factors would be, at best, a 
guess. 

Three surprising features were observed in simulation 
deposit patterns. ‘The first, the apparent rotation of the 
(100) pattern at 1 keV, will be further discussed. The 


rotation of the (100) pattern for low energy sputtering has 


not been obServed. ‘The results of the simulation indicate 
that such a rotation is possible. This conclusion is based on 
empirical rather than theoretical considerations. ‘The hypocycloid 


outline found at 3 keV in the simulation appears to be formed 
by the same atoms which form (211) spots at 1 keV. 
Furthermore, no explanations have been previously proposed as 
to why the outline is a hypocycloid in experimental patterns. 
One would expect that a circular or perhaps a non-distinct 
outline would be observed when the (100) surface is sputtered. 
Accordingly, it is proposed that the hypocycloid outline is a 
result of the inability of the (100) surface to completely 


focus atoms sputtered by a Squeeze mechanism for bombardment 
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energies at which the energy propagation 1s not contined in the 
surface. 

The second feature was the unique form ot the bonds which 
connected the spots in the simulation pattern for (111) 
copper sputtered at 5 keV. ‘The appearance of the similar 
outline, a spherical triangle, was obServed in a pattern of (111) 
copper sputtered with 1.5 keV krypton by Yurasova and 
Bukhanov(34), The third feature was the hexagon spot pattern 


(10) Mong 


at 2 keV, a pattern Seen by Anderson and Wehner 
(111) copper sputtered with 400 ev mercury. The similarities 
observed in the simulation patterns with theSe anomalous 
features of the experimental] patterns Suggested that common 
factor other than the (111) copper surface might be present. 
The momentum ratios of 400 ev mercury=2 keV argon and 
15keV krypton=-5 keV argon are 1.0 and 0.8 respectively. 
This is not considered to definitely establish a sputtering 
correlation based only on ion momentum. ‘The momentum ratiow 
for 1.5 keV krypton and 4 and 3 keV argon are 0.9 and 1.0 
respectively. ldeally, the spherical triangle would be observed 
in the 3 keV (111) simulation pattern. Additionally, when 
(111) copper is sputtered at high temperatures, the triangula? 
outline has been observed to become more pronounced (34) , 
The absence of such @€ momentum scaling ettfect coujid be 

easily shown by experimentally sputtering (11') copper with 

liom Momentum compatable to that of 400 ev mercury. 
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The general consistency of the numerical values of 
sputtering ratios, and energy distributions of sputtered atoms 
and likeness of the patterns with experimental data is a 
strong argument for the validity of a computer simulation of 
sputtering. It is remarkable that the consistency is as 
good as it is. The model uses only a repulsive potential, the 
crystal size used in the simulation is an infinitesimal portion 
of the smallest laboratory specimens, and the potential form 
and parameters are comparatively crude. It is concluded 
that the results obtained substantiate the concept of 
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APPENDIX A 
The Beam Model 
ee apace areas and Pimpact points. 

Fiach crystal surface contains an intrinsic, plane geometric 
shape; hexagon for (111) surface, square for (100) surface, 
and rectangle for (110) surface. A volume element of the 
(hkl) surface is defined by the area of this intrinsic shape ; 
and a depth of some number of (hkl) planes. This volume is 
chosen so that by translation along axes of a Cartesian 
coordinate system the entire crystal may be generated. The 
smallest intrinsic area which can be chosen for each surface 
and still satisfy the translation requirement is shown in 
figure 72. A finite number of points is symmetrically distributed 
within each of these areas to represent the infinite set of 
possible points of impact for an incident ion. ‘These smallest 
aqetowarce i urener divided Into representarive imMpace areas 
which are Seen from figure 73 to be degenerate under appro- 
priate rotations and/or inversions of the coordinate axes. 
=mce ere seceors are degenerate, the ser or imac cmn oto 
in each Sector is also degenerate; only the points contained | 
in one impact area need be used to represent bombardment of 
the entire area. The coordinate axes rotation and inversion 
schemes are discussed in Appendix C (Deposit Pattern 
Produc zion). 

The independence of the impact point set used, with 
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respect to number and Kinetic energy of atoms sputtered, was 
tested using two sets of points. No dependence was found. 
Results for the (111) surface using a 3 keV ion are shown in 
figure 74 as an example. Set 1 points are those shown in 
figure 73, set 2 (not shown) is a set of eleven points which 
are located between the points of set 1. 

B. Neutralization of Beam ions. 

The assumption that argon ions are neutralized prior to 
impact on copper may be inferred from results of a theoretica) 
study of secondary electron emission by Harrison et. al, e0), 
Consideration of atomeatom rather than foneatum interactions 
gave results in reasonable agreement with experimental data. 

A supporting argument for neutralization is based on 
Hagstrum's theory of Auger ejection of placenora ale ‘The 
probability of an ion Deing, neutralized in dx at x is: 

Py. (x,v) #a exp) ~~ exp [a(x )| ma(x=x_ ) 
where Xm = (1/a) Qn (A/av) is the value of x White dw 1s oa 
maximum. ‘I'he parameters A and a occur in the wnesiaee ion 
rate function, and v is the velocity of an ion starting at x = Oo. 
Hagstrum used tungsten as an example and obtained the vali 
of A by empirical means, the value of a from published data. 
Rather than attempt an exact proot of neutralization tor argon 
on copper, Hagstrum's results tor argon on tungsten have been 
used to give an order ot magnitude, at: worst, approximation. 

S 


Accordingly, 1.632 xe, a*.oae ror] to 10 keV argbi ions on 
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copper; this is a reasonable distance from the crystal surface 
and one may assume that neutralization occurs. _ 
Ge lon WMeflection by Surface Potential. 

As a first approximation, it is assumed that a singly 
positive-charged argon ion sees the crystal as a single, fixed 
scattering center having Zot positive charge. The well known 
equation for a central force-induced hyperbolic orbit is used 
in conjunction with figure 75 to determine the deflection. 

r(Q@) = (L*£/mk) Pct + (14 2EL2/mK2)1/* cos( @- Oe | 

tL = ion angular momentum m = ion mass 
E 


© 


, = angle of closest approach 


I] 
H 


ion energy ine Q, 2/4 WE. 

The angular momentum L is determined at ro, O97 
to be L = mvb where mv</2 = & and DMIs@Ehe immace parameter. 
The equation for r(Q) is put in a more Gonveniente formubes 
making the substitution we = (mvb)* = 2mEb*~ to give 

r (QO) =b (2Bb/)/f —1 + E (28b/k)® | cos(O-@, ) 
The angle of closest approach Q, is determined as a function 
Of the product Eb by Pequiring the denominator to vaio 
O=77 . Once O, (as “ar fune tion =o Eb) is known, the 
deflection ratio alo is calculated at the crystal surface. 
The deflection calculated is for an unneutralized ion, thus 
greater than that for an ion which 1s neutralized at Some 
distance in front of the surface. Plots of r/b at the surface 


and O as functions of the ion kinetic energy-impact parameter 
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product are contained in tigure 76. The numerical values are 
for an Ar’ — Cu system. 

A more sophisticated approximation, scattering from a 
fixed dipole of same charge sign, requires a messy integration. 
Rather than follow this line, one may expect from the nature 
of the problem that the path of the approaching ion will undergo 
some oscillatory motion or perhaps corkscrew motion as the ion 
is influenced by the surtace potential. In either case one 
would expect that the net acceleration of the ion parallel to 
the crystal surface would be no greater than that due to a 
single fixed scattering center. With these considerations, 
figure 76 is used to determine the percentage of ions which 
will be appreciably detlected. As an example of an appreciable 
ion deflection, consider r/b = 2.0; from figure 76, Eb = 3.5 x 

O 
iene keV-A, An energy range of 1-10 keV Yor argon ions, 
that used in this study, corresponds to & range of impact 
O 

parameters of 0.035 = 0.0035 A. ‘The fraction of surface 
area, and theretore traction of fons which will be appreciably 
deflected is (0.035)*/(1.26)" = 7.7 x LUT? for 1 keV ions ana 

) 
7 es 1078 for. iv WeV Jjionsealchdee A is the copper 4&tam 
effective radius in the crystal). ‘This negligibly small traction 
of ions cannot influence the macro aspects of the sputtering 
model. 
D. Equilibrium State During linpact. 


The use of a single atom approximation to the beam with 
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the crystal model described in the main text requires that the 
crystal region be in an equilibrium state at the time of an 
ion impact. It is not necessary that this equilibrium state 
be identical to previous equilibrium states since the eueioiey 
of surface configurations available in the model provides for 
random surface conditions. It is necessary that the time 
required for the crystal region to return to an equilibrium 
state be small with respect to the arrival-time intervals of 
the ions. Satisfaction of this condition is determined by 
comparing an experimental beam flux over the area of the 
crystal face of the model used to the time required for 
completion of all energetic collisions in the model. 


a were used 


Beam intensities of the order of 1004 amps/cm 
by Magnuson and Carlston!12) , If a beam of this intensity 
is incident on a crystal surface area of less than 1032A* such 
as in the model, the ion flux over this area is less than 100 
ions/sec, an ion arrival-time interval of suas seconds; all 
energetic collisions in the model are completed within about 


1 seconds. The relaxation time of the region is thus much 


10 
smaller than the ion arrival-time intervals and the condition, 


that the crystal region be in an equilibrium state at the time 


of impact, is satisfied. 


re 


AS aaeN Dix BB 
Positioning the lon 

The ion is positioned tangent to the first target atom it 
will strike. Since this is a dynamic rather than static process 
it is not necessary that this be a stable position on the crystal 
surface. Figure 77 shows the ion at its arbitrary initial 
position and calculated final position, both with respect to the 
impact point and target atom. ‘The initial position is a small 
distance in front of the surface, beyond the eroded potential 
range of the crystal atoms. ‘T’his position is described by a 


vector ry originating at the impact point and having direction 


parallel but opposite to the ions velocity vector. ‘The vector 


ro, from the impact point to the target atom, is known 


Sijge Uae Gane acon secoordinatesware known; the vector 1. 


is known, and it is desired that vector r have a magnitude 


cout 


equal to the distance between centers of two tangent atoms, 


ar,- time law of cosines 1s used to compute the magnitude of 


oS ae 


a which lies along rae Accordingly, the following calculations 


—Saa 
Cee an 


are made: 


cosa@= Yr .,Y 
1 = 
ra} ral 
(r3) - (er, cose ) (r,) 4 oss si (er) ene 
Yr eo . » Q : ee . «2 a 6 h49 \& ° : 
(v3) r, Cos E cos & - (2, = (er) ) ie 


The positive square root solution is chosen to give the tangent 


on the outside hemisphere of the ltlarget atom. 


72 


The ion positioning is accomplished automatically for each 
run by subroutine START, which is contained in the computer 


program for the sputtering simulation. 


i. 


Pee NAD TK aC 

Production and Analysis of Sputtering Deposit Patterns 
2, Se gmodmction. 

A sputtering deposit pattern represents the intersection 
points of atoms! velocity vectors with the surface of a 
collector plate. The sputtered atoms from the simulation 
are collected on a flat plate by determining thesSe points of 
intersection. Bach atom which exits through the crystal 
surface is initially assumed to have been sputtered and a data 
card has been prepared for each onewes Data of maxticular 
interest are each velocity component magnitude and the Kinetic 
energy perpendicular to the surface; data of secondary 
interest are atom number, impact point used and ion kinetic 
energy. (The use of these last data will be discussed shortly.) 

It is recognized that each sputtered atom has lost some 
energy to overcome the surface binding energy. An assumed 
value of binding energy is subtracted from the perpendicular 
kinetic energy and a new perpendicular velocity component is 
calculated. Parallel velocity components are normalized to the 
new perpendicular component to give a two dimensional coordi:- 
ate point. ‘This point is the intersection of the atom's 
velocity vector with an imaginary collector plate placed at unit 
distance from the target surface. Hach point is then rotated 
and/or mirrored about the coordinate axes to give the inter- 


section of a velocity vector which would have resulted for an 


Th 


ion impacting at the corresponding impact point in each of 

the other impact areas. The impact areas are shown in 
figure 73 and the coordinate point rotation and mirroring values 
are listed o program DATASORT which is used to generate 
Ele Ppolntce == “ach point is plotted Using program DATA LOT 
with a CDC1i60A computer and a CalComp plotter to give a 
point plot. ‘The dimensions of the plots are in target—to- 
collector or T-C units since these points have been normalized 
to unit target-to-collector distance. The scale which has 
been used permits plotting deposit points of atoms which have 
been sputtered within an escape cone of about 63 degrees 

(57 degree cone shown in figures). This has been found 
Satisfactory to contain all pattern features of interest. 

The point plots show only point patterns and therefore 
do not accurately simulate experimental patterns which are 
eames patterns. ‘lhe conversion Of a pelle plot ued 
smooth area density pattern is made by photographing the point 
plot with the camera defocused such that no single point is 
distinguishable but high and low density areas are prominent. 
Developing and printing is controlled to bring out the high 
density areas while maintaining the haze background. Loss of 
intensity in Some spot regions cannot be avoided such as Seen 
in figures 6b and 12b. The process must be adjusted for 
each pattern; the Sequence of photography, developing, and 


the printing of the positive image is highly dependent on the 


CS, 


ability of the photographer. Concurrently, the quality of 
reproduction of these patterns in printed torr depends on thie 
plate preparation, the printing press, and the paper used. 
B. Analysis. 

A 30 x 30 square grid is placed over the central 3.0 x 3.0 
unit square of the raw pattern by program DATAGRID. 
The identification number of an atom and the impact point and 
impact area ot the run in whieh that atom was sputtered are 
recorded for each atom in the grid square. The total number 
of atoms, the total energy, and the average energy per atom 
for each grid square are printed in separate Square arrays. 
The individual grid square data provides for correlation 
between an atom's crystal location and its deposit point in 
the pattern; the square arrays of the number of atoms and 
the total] and average energy densities may be compared 
directly with the point plot or smooth pattern for analysis of 


pattern features. 
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Secrion »..: 


Section. —: 


Section 3: 


APPENDIX © 
Memory block allocations and functions are 
established. All storage cellS are zeroed. 
Input and output formats are specified. 
Constants are set, and the target material, 
ion Species, and crystal face to be sputtered 
are read in. 
The ion energy, impact area, and impact 
point are read in. Constants peculiar to the 
run are eStablished. The appropriate sub- 
routine is called to build the crystal, and 
crystal boundaries are calculated. Subroutine 
START is called to position the ion. The 
initial value of At is calculated. Initial 
coordinates of all atoms are assigned, and 
their velocity components are zeroed (except 
for thevion). Lwelinitial’ coordimaees ct the 
crystal atoms are printed. 
Forces are calculated by calling subroutine 
STEP. Atoms are moved to their intermediate 
positions in the two step cycle. Subroutire 
STEP is called again, and atoms are moved 
to their final positions. Final velocities are 
then computed. Force components are 
zeroed in preparation for the next timestep. 


ew 


The maximum kinetic energy is determined for 
the calculation of At. 

Section 4: The time remaining until cutoff is determined. 
If there is insufficient time to complete 
another timestep, terminal data is printed. 
Potential energy is calculated and summed 
with kinetic energy to give the total energy 
for energy balance check (manual check). The 
data for atoms having potential energy greater 
than the thermal energy are printed. 

Section 5: Atoms which have kinetic energy but not 
potential energy are assumed to be free of 
the crystal. They are assigned LCUT=1. 
The surface through which an atom exited is 
determined and a code assigned. A maximum 
potential energy is found (for &t calculation) 
among atoms Which do not have LCUT=1. If 
the maximum potential energy is less than a 
minimum value the terminating process begins. 
If this energy is greater than the minimum 
value, another timestep begins. 

Section 6: Pertinent data for all atoms is printed. Atoms 
which have exited through the front of the 
crystal or will exit through the front are 
assumed to be sputtered. Data for these 


10). 


atoms areprinted and a data card punched 
for each atom. A new data card is then 


read into initiate anette rman. 
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COX! 
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ez 
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CVD 
oy 
CN 
ey 
weox 
DeoOy 
Deo zs, 
DEE 


DIST 


ZS a as, 


CGleoeeamy forek CCSPrU'T 


Parameter for target force function correction. 
Pieating spomiemetorm of [xX > LY ,i]. 


Cosine of the angle between vectors: Ri, Re. 
ALFA squared. 
Reciprocal of magnitude of atom velocity. 


Energy which an atom within the crystal at shut-. 
down must have to be considered sputtered. 


Variable representing primary material. 


Unscaled x, z coordinates of the impact point 
Frictional force multiplier. (See CVS) 


Direction cosines of primary velocity vector. 


Negative values of COX, COY, CO2&. 


A constant. 
Converts meters to angstrom units. 
Converts electron volts to joules. 


Converts atomic mass units to kilograms. 


Converts lattice units to angstrom units. 
Direction cosines of sputtered atom velocity vector. 


Distance difference between nearest neighbor 
distance and actual atom differences. 
Distance between any two atoms. 
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EMAX 
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Vp os) 


POR 


xX, V, 2 COMpoenenes or 2h. 2. 


Length of timestep in seconds. 


Number of lattice units most energetic atom 
may move in one timestep. 


DT/CVD - a ratio used to avoid repeated division. 


DT/PTMAS - a ratio used to avoid repeated 
division. 


DT/PGMAS - a ratio used to avoid repeated 
division. 

x, y, z distances atom has moved from initial 
position. 

A lower limit on an atom's potential energy. If 
energy is less than or equal to ECUT the program 
shuts down. 

7S CUCOll veneers yy. 

The maximum energy encountered in any cycle. 
Primary energy in electron volts. 


Primary energy in kilo-electron volts. 


Potential function parameters. 


The component force increment on an atom. 


The minimum distance the primary is positioned 
in front of the first xz plane at start time. 


A. small number used in checking potential energy 
ZCrO =peini. 


FORCE/DIST = a ratio used to avoid repeated 
division. 


Numerical valve Cf “tne £orce fumetsion with a 
variable parameter. 


LOS 


FORE 
FPTC 
FPTE 
FRC 
EX 
FY 
F2 


EF XA - 


GMAS 


EWiDar© D 


PDTOM 


EO TOMB 


RentAS 


reer SS 


eetsT 
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iS 
cs . 
iL 
liga 


IN 


LIND eho a4 


INOW 


Target atom force function. 

The corrective force value at ROE. 

‘The corrective force function. 

Numerical value of the target force function at 
Ree... 


xX, VY, 2 COmponents of total force on an aisom. 


Force function parameter. 
Target atom mass (in a.m.u.) 


1/2 DTOD - a ratio used to avoid repeated 
division. 


1/2 DTOM - a ratio used to avoid repeated 
division. 


1/2 DTOMB - a ratio used to avoid repeated 
division. 


1/2 GMAS - a ratio used to avoid repeated 
division. 


1/2 TMAS - a ratio used to avoid repeated 
division. 


Used to provide output prior to time limit 
shutdown. 


Dummy variable. 


Alpha-numeric arrays for titling. 


Odd-even integer used to determine atom site 
establishment. 


Integer (0 or 1) used in determining dynamics 
cycle step. 


Time program has been running in seconds. 
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die 


Ji Real 


3,6 
stg 
Le, 
ae? 
HSS alee 
Beale 


‘relat Wal 


Smal 


eS 
earl 


re 7 TF 


Gal 


eC UT 


Subscript value of atom. Wsed@tin subreurige 
Sie 


Time left prior to time [limic. 


Fixed point valiies 01 "=P 20, 9] he. 


Subscript value of most energetic atom. 
Unscaled fixed point x coordinate used in lattice 
generation. Also a dummy variable in function 
ea 

Odd-even integer used to determine atom site 


establishment. 


NUM DerMer cy" Vy mee plates Gleascig Sital . 


Grystal dimensions in=say,= 2 - 

Cutoff variable based on total potential energy 
Of “Crystal: 

Unscaled y coordinate used in crystal generation. 


Variables used to establish atom sites. 


Identifies exit point of atom. 
Unscaled 2 coordimaeemiece to establish atemiecite. 


Used to identify atoms which are not included in 
calculations. 


The highest numbered atom in the crystal. 
Sum of the Miller index integers. 


Used to identify type of surface, i.e., regular, 
Stub, vacancy. 


An integer used to begin atom numbering. 
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Subscript of the atom to be removed for 
Vacancy ssumrrace. 


Data SoOuiepiie. Incremecin. 

Page numbering variable. 

Initial print statement cycle. 

Cutoff variable based on too long a timestep. 
Timestep. 

Timestep limit before shutdown. 

Same as AC except applicable to primary. 


Primary force function parameters. 


Pepeconstant’ = 0,53 

Pritaryeflorece. functions 

PEFPTE evaluated at ROE. 

Frimemry corrective force Tunction. 
BORE -evaluated atreROkE. 

Peunary force funetionsparametGer: 
Primary mass in kilograms. 

Kinetic energy of an atom. 

Y component of kinetic energy of an atom. 
Crystal plane (alphanumeric variable). 
Same as PLA. 

Impact point (alphanumeric variable). 
Potential energy between two atoms. 
Target potential function. 

Potential energy of an atom. 
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Beimary potential function: 

PRPOTE evaluated ae, mee. 

Chemical symbol for primary maeerial. 

POTE evaluatedwabe. ei. 

Total energy of an atom (potential kinetic). 
Target mass in kilograms. 


A small number used in checking kinetic energy 
ZeELOLpOlIe. 


Cutoff variable checked against total potential 
energy. 


Unecaled x, 2 COordinaveswor=riipact area 


reference point 
Nearest neighbor distance. 


ROE - DY) (one timeseepeetouance less than 
Miatie Cisteamiee |r, 


ROE squared. 


X, YY, Z2 €eccrdigaeescworeatom at any time, 


X, Y, Z coordinates of crystal boundaries other 
than zero. 


xX, ¥Y, 2 COordinates ov “atoms initial position: 


x, ¥, 2eCCCeiaeescee LCemporary position of 
atom during foree cyele. 


my V> 2 ceoemeireeeset linpact point. 


VeCCOP fol iret etoolteeLo initia) primary, Besition. 


ils 


Foe 
ey, 
ees 
Pere 


Fe a 


Ee 


3 


aa 
SCX 
= 
Soe, 
SLOW 


=e 
ae 


SS 2, 


SOS 


IANS 


a G Er 


lag al 


“10 a 


ae be IM 


cae 


eM E; 


ae PIO 


TMAS 


Xo Vig 2 COOrdinates of INibialermimary position. 


Scalar product Of VeCtors =. 


Magnitude of vector from impact point to first 
atom hit by primary. 


Re squared. 


Magnitude of vector from impact point of 
primary start position. 


i eeeOT . 
me y, 24 coordinate scale factors. 


Cutoff variable checked against a long DT. 

4 
x, z distance from impact area reference point 
GOmlimpacs polmt. 


A z scale factor used for (111) plane Mattice 
generation. 


Plane (alphanumeric variable). 
Chemical symbol for target material. 
Target material (alphanumeric variable). 


Total energy of crystal atoms (kinetic + 
potential). 


Temperature of lattice in degrees Kelvin. 
Thermal energy of atom. 

Computer time program has been running. 
Elapsed problem time. 

A. function to convert seconds to minmweecs 


Target atom mass in kilograms. 
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Total kinetic energy of crystal atoms. 
Total potential energy of crystal atoms. 
Magnitude of primary velocity vector. 


Storage variable for velocity components. 


X,.¥, 2 Compemencs of atom velocity. 


Unscaled x coordinate used in crystal generation. 
EFleating point form of 1D weet. in function 2iie 
Unscaled y coordinate used in crystal generation. 


Unscaled z coordinate used in crystal generation. 
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F(r) 


Erosion of the force function 

atr= 2X6 results ina loss 

of energy in the model. This 

is corrected in the STEP subroutine. 
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Figure 15b 
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Figure 16b 
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19. ode 8 : 

e sputtering process has been investigated by simulating the 
sputtering of single-crystal copper with 1-7 keV argon. A digital 
computer was used to build the crystal, bombard it, and move 
crysta] atoms. Four mechanisms were observed which cause surface 
atoms to sputter. An atom is sputtered when (1) it is squeezed 
out of the surface, (2) it ig scooped out when another atom strikes 
its inner hemisphere, (3) it is ejected when an atom passes behind it, 
and (4) it is knocked out by a second layer atom which is moving 
outward. Wearly all sputtered atoms were surface atoms. Second 
ang third layer atoms were sputtered only for ion energies greater 
than § keV. They were sputtered by mechanisms similar to the 
surface atom mechanisms. "Silsbee chains" were observed to be 

| directed into the crystal, and momentum focusing was observed to 

| cause sputtering only when it occurred in close packed, surface rows. 
Qutward directe@ chains were not observed. Sputtering deposit | 
patterns, sputtering ratios, and sputtered atom energy distributions 
were obtained for (100), (110), and (111) surfaces. All data 
compared favorably with experimental data. 
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